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Abstract 

"-H , 

p^ ' We consider the modes of the electric field of a cavity where there is 

\^ • 
j^ ' an embedded polarized dielectric film. The model consists in the Maxwell 

equations coupled to a Duffing oscillator for the film which we assume 

infinitely thin. We derive the normal modes of the system and show that 

they are orthogonal with a special scalar product which we introduce. 

These modes are well suited to describe the system even for a film of 

finite thickness. By acting on the film we demonstrate switching from one 

cavity mode to another. Since the system is linear, little energy is needed 

for this conversion. Moreover the amplitude equations describe very well 



this complex system under different perturbations (damping, forcing and 
nonlinearity) with very few modes. These results are very general and 
can be applied to different situations like for an atom in a cavity or a 
Josephson junction in a capacitor and this could be very useful for many 
nano-physics applications. 

1 Introduction 

There are many reasons to couple an oscillator to a cavity. One example is a laser 
built using a Fabry-Perot resonator enclosing an active medium which can be 
modeled as a two level atom [U [5] . The cavity can also be used to synchonize 
oscillators [3] as for an array of Josephson junctions. For window Josephson 
junctions, used as microwave generators, the Josephson junction collects all 
the energy in one of the cavity modes pj. The cavity can also be used as a 
thermostat to cool the oscillator as described in [S] where an optical cavity 
is used to cool an atom. In another example coupling an atom or quantum 
oscillator to a resonator can significantly change its transport properties |6J . 

In all the systems described above we have a localized oscillator coupled 
to a resonator. In addition the size of the oscillator can often be neglected. 
This situation can be represented by a thin film model [7l [8] . Such a thin film 
model was generalized by taking into account the local field effects (dipole-dipole 
interaction) [Hllinilll]- Intrinsic optical bistability is the main result of this gen- 
eralization. Thin films containing three-level atoms [12_, two-photon resonant 
atoms [131 m] , inhomogeneously broadened two- level |TS] atoms and two-level 
atoms with permanent dipole [16] represent the different generalizations of the 
model. The coherent responses of resonant atoms of a thin film to short optical 
pulse excitation were considered in [15] . It was shown that for a certain intensity 
the incident pulse generates sharp spikes in the transmitted radiation. Photon 
echo in the form of multiple responses to a double or triple pulse excitation 
was predicted also in this paper. The coherent reflection from a thin film as 
superradiation was studied in [lH [ITl [19] . 

Recently we used the thin film model to describe switching phenomena in 
ferroelectrics 20 . The behavior of the new artificial materials - metamaterials 
- under action of electromagnetic pulses could be also described by this model 
|21| . As we see the thin film model is extremely fruitful, the investigation of the 



behavior of the thin fihn embedded inside the cavity is a very attractive prob- 
lem. When the model of thin film is explored, the problem of the matter-field 
interaction reduces to the one of an electromagnetic cavity with an embedded 
(linear or nonlinear) oscillator. Frequently the nonlinear problem is analyzed 
using coupled mode equations. In this approach the relevant values are the am- 
plitudes of each of the linear modes. The nonlinear partial differential equation 
of motion is reduced to a system of coupled ordinary differential equations for 
the amplitudes of the linear modes. A suitable choice of the linear modes allows 
to get an effective description of the original problem. 

We consider here this simplest and most general situation, first for small 
energy for which the problem is linear. In this case, the system film/ cavity has 
normal modes which we characterize. These describe correctly the evolution of 
the system as opposed to standard Fourier modes or even other simpler normal 
modes. These adapted normal modes are orthogonal with respect to an adapted 
scalar product. Using these modes we can define simply the state of the system 
and obtain mode conversion by acting on the film sub-system. In this paper 
we propose the mode conversion mechanism for a linear system which can be 
used for different applications. In the nonlinear case, for medium amplitudes, 
we show that a few modes are sufficient to describe the evolution of the system. 
After introducing the model in section 2, we consider the linear limit in section 
3 and derive the normal modes of the system. The special scalar product is 
derived in section 4. In sections 5 and 6 we use the normal modes to define 
the state of the (linear) system and show mode conversion when driving and 
damping the film. We also describe the general nonlinear case and we conclude 
in section 7. 

2 The thin film model 

We consider a one dimensional model of the electromagnetic radiation interact- 
ing with a polarized dielectric film inside a cavity. The film is placed at the 
distance Xa inside the cavity having the length /. The Lagrangian density for 
the electromagnetic field, the film medium and their coupling is the following 

«t "? C/ N f 1^ <f 9^ \ ,-,^ 

C = — - — + 5[x-Xa)\—-m—- — -aqat\. (1) 



Here a is the analog of vector potential and q is the medium polarization, a 
is a coupling constant. The last term in Lagrangian describes the coupling 
between a and q. The dielectric medium can be ferroelectric {m = — 1) with 
two polarizations or paraelectric m ^ 1. The Hamiltonian of the system is 



which gives 



dC dC „ ,„, 

n^at--+qt—-C, 2 

oat oqt 



H = | + | + 5(x-.j(| + m^ + ^^), (3) 



The variation of the action functional yields the Euler-Lagrange equations for 

a and q 

dC d dC d dC , , 

da dt dat dx dax 

dC d dC d dC . , 

dq dt dqt dx dq^ 

which reduce to 

att-a^x = aS{x-Xa)qt (6) 

qtt +mq + q^ = - aat (7) 

The equations for the electric field e = —at and medium variable can then be 
obtained 

ett-exx = -aS{x~ Xa)qtt, (8) 

qtt + mq + q^ = ae{xa), (9) 

where the coupling between the fields e and q only occurs in the medium at 

X — X(^ . 

In a recent article [5D], we considered with this model the interaction of a 
thin dielectric film with an electromagnetic pulse. We studied both the case 
of a ferroelectric and paraelectric film. For the ferroelectric film we showed 
that the polarization can be switched by an incoming pulse and studied this 
phenomenon. Here we will assume that the film is embedded in a cavity and we 
will study how cavity modes can be controlled by the film. Specifically we will 
assume Dirichlet boundary conditions for the field. 

For small amplitudes of the field, it is natural to neglect the nonlinear re- 
sponse of the film. Note however that the ferroelectric film and paraelectric film 



have different naturai frequencies corresponding to different stationary points. 
For the paraelectric case, there is only one fixed point q = while for the fer- 
roelectric case there are three fixed points, the unstable one q = and the two 
stable ones q = ±1 corresponding to two opposite signs of the polarization. It 
is then natural to introduce the natural frequency of the oscillator luq — m for 
the paraelectric case and ujq — m + Sgp for the ferroelectric case. We therefore 
consider below the general linear problem of a harmonic oscillator of frequency 
Wo embedded in a cavity. 

3 The linear limit: normal modes 

The linear problem is 

ett-exx = -aS {x - Xa) qtt, (10) 

qtt+ujlq = ae{xa). (11) 

Note that we have a Dirac delta function in the first equation so that the solution 
will not have a second derivative at x — Xa- In this case, one can write the 
solution using standard sine Fourier modes. However these are not adapted to 
describe the evolution because the projection of the operator gives wrong results 
[22]. Then we need to define new normal modes. For this one first separates 
time and space and one looks for solutions in the form 

e(x, t) = E{x)e''^\ q{x, t) = Q{x)e"^' 

so that the system pUj) becomes 

E"{x)+u;^E{x) = -auj^QS{x~Xa), Q^^P^, (12) 

LOq — to 

As expected from the general theory of linear operators [,23^ the system will 
exhibit eigenfrequencies and eigenmodes (normal modes). Combining these two 
equations, we obtain the final boundary value problem for E 

E"{x) + c.^ ("l + "'^^"°^ ) E{x) ^ 0, (13) 

with the boundary conditions E{Q) = E{1) — 0. 

To obtain the solution, notice that except for x = Xa 

E"{x)+LJ^E{x) = 0. 



Using this remark and the boundary conditions we get the left and right sohi- 

tions 

, , I Asinwa:, x <Xa-, , ^ 

E{x) = (14) 

(^ BiiinUJ[l — X), X > Xa, 

where A and B are constants. To connect the left and right solutions we use 
the continuity of E{x) as well as of e at x = Xq. The second relation needed 
is the jump condition for E' obtained by integrating (|13p over a small interval 
centered on Xa ■ When the size of the interval goes to zero we get 

{E']S = ^-f^E{x^). (15) 

At X — Xa the continuity of E and jump condition (jlSp give the following 
relations 

AsiiiujXa — B si'iiuj{l — Xa) = 0, (16) 

2 

A(—^ TT smujXa — cosujXa) — B cosujU — Xa) — 0. (17) 

Wq — UJ'^ 

For this homogeneous system to have a non zero solution, the determinant must 
be zero and this gives the dispersion relation 

smu!l = a'^^ ^ sinijja;a sina;(Z — Xa), (18) 

UJq — LO'^ 

which determines the allowed frequencies lo of the system. 

As can be expected from the general theory [53], we have an infinite count- 
able set of allowed frequencies. Note that in the absence of coupling to the film 
a = 0, we get the standard Fourier modes cij„ = wk/I. For small a the shift 
in frequency is small because the right hand side is proportional to a^ . These 
eigenfrequencies can be computed using bisection for example. 

To summarize, the eigenvalues and eigenvectors (up to a multiplicative con- 
stant) of the boundary value problem p2p are 

where 



sm LOiX 
sin uj 



jf^^SUYLOi{l- X), X>Xa, 



and uJi satisfies ((TH]). 

Note that an equation similar to P^ was obtained in the theory of a ID 
waveguide with a perfect mirror at one end and a two-level atom at the other 
end [24]. Contrary to our case, this is not an eigenvalue problem because the 
system is open on one side. 

From another point of view, the system cavity/film (fTU|) was considered for 
Wo = by Bocchieri et al [U in the context of statistical mechanics. Their 
main result was that there was always energy exchange between the film and 
the cavity so that equipartition could not be reached. Indeed, this can be seen 
by examining the normal modes (j20p which couple Ei and Qi. Only for special 
symetries (film at the center of the cavity . . . ) and special frequencies do we 
get Qi = 0. 

3.1 Influence of the film parameters on the dispersion re- 
lation 

We now study the dispersion relation in more detail. In Fig. [1] we plot the 
solutions of (fTSl) as a function of the film position for t^o = 1 and a = 1. 
Notice how the systems generates two eigenvalues in place of the frequency 
loq = I. For large frequencies u '^ loq we recover the standard Fourier cavity 
modes ujn = nir/l. The eigenmodes for the particular case of a centered film 
shown in Fig. [1] contain the even Fourier modes. These correspond to Q^ = 
because their derivative is continuous at Xa- This will have important practical 
consequences. 

Because of the film, the eigenfrequencies of the system film/cavity differ 
from the usual Fourier cavity modes uj„ = nn/l. They are shifted if tOn 7^ ^o 
and disappear if uJn = loq. Let us compute this shift in the limit of small a 
using perturbation theory. To simplify the analysis, we assume Z = tt so that 
UJn — nil /I — n. We search the frequency uj using the expansion 

u! = n + a LJi + . . 

with 071 ^ 1 and n are the usual sine Fourier modes of the cavity. We have the 
following relations 

smujl^ (-I)"a2wi7r + 0(a^) 

sinutXa = sinna;^ + a uJiXa cos nXa + 0{a ) 
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Figure 1: Eigenvalues (zeros of the dispersion relation ^TE\i ) as a function of 
the position of the film Xa- The parameters are I — tt, a — 1 and ujq = 1 



sina;(/ — Xa) = sinn{l — Xa) + a llJi{1 — Xa) cosn{l — Xa) + 0{a ) 
Plugging these relations into (fT8|l we have, up to 0{a'^) 



(-l)"^i 



7r(w2 _ „2) 



sin nXa sin n(l — Xa)- 



Assuming that a^ < 1 we get the simplified expression 

uj = n — a (— 1)" o ^ sinnxa smn(l — Xa) + 0(a ). 



(21) 



(22) 



Due to the presence of the film the cavity modes close to the film mode are blue 
shifted if the frequency of cavity mode is above the oscillator eigenfrequency 
and red shifted for lower cavity mode frequencies. 

When ljq = tiq where no is an integer, the oscillator frequency coincides with 
the cavity mode. In this case the eigenfrequency of the combined system splits 
away from hq. Again this can be calculated for small a by assuming 



UJ = UjQ +UJI, 



where \uji\ <C uJo- 



Plugging this expression into the dispersion relation and collecting the terms, 
we obtain the second degree equation 

Alui"^ + Blui + C = 0, 
where 

A = 47rn + a^((7r- 2xa)sm2nxa - {cos2nXa - l)"/2), (23) 

B — a^(n(7r — 2xa) sm2nxa + cos2na;(j — 1), (24) 

C = a^n(cos2na;a - 1). (25) 

The two branches of the resonant frequency are then given by 



-B ± VS2 - AAC 

^^ = 2A • 

As an example, consider the case of Fig [T] corresponding to / = tt and a film 
placed in the center of the cavity Xa = 0.51. Then cosa;(/ — 2a;a) = 1 and 
cosul^ (-l)"o(l-wJ/2) 

Then the splitting is given by 

AttuoujI ^ a^m + uji){l ~ ujI/2){1 ~ (-1)"°) (26) 

so that wi = for even resonant frequency uq and there is no shift from the 
resonance in this case. For the odd rip we obtain the quadratic equation for the 
frequency detunings uji from the resonance: 

"' =0 (27) 



no(27r + a2/2) 2T: + a^/2 



4no2(27r + aV2) 
2n(27r + a72) ' 



^1 = o^.o, , ^2/o^ I 1 ± \/ 1 + .2 — 1 (28) 



When a increases, the influence of the film grows and the eigenfrequencies 
become very different from the Fourier cavity modes. In fact when a ^ 1, the 
dominating term in the dispersion relation is the right hand side and we obtain 

sinwxa = 0, or siiLuj{l — Xa) — 0- 

so that 

mr mr . 

LOn = — , or w„ = — . (29) 
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Figure 2: Plot of the first fO eigenfrequencies as a function of the couphng 
parameter a for an off-centered film Xa ^ 1/4 = ir/A whose frequency is loq = 3. 

This corresponds to oscillations in the left or right cavities. Fig. [2] shows the 
first 10 eigenfrequencies as a function of the coupling parameter a for a cavity of 
length / ~ n and an off-centered film Xa — //4 whose frequency loq = 3. Notice 
the splitting for uj — 3. As a increases, the eigenfrequencies tend to the ones 
given by (|29|) . ie w„ = in or w„ = 4n/3. 

4 Orthogonality of the normal modes 

Using the vector notation V defined above, the original linear system ^TU\i can 
be formally written as 



where the operator L is 



a'^5{x — Xa) 



aujQS{x - Xa) 
2 



,0 Oy \—a J d{x — Xa) ujQ 

The eigenfrequencies and eigenvectors uji, V^ are such that 



(30) 



(31) 



LY.^LufY,. 



(32) 



The boundary value problem (|13p is not a standard Sturm-Liouville problem 
because the potential depends on the eigenvalue. Therefore one needs to define 
a particular scalar product so that the normal modes V^ defined previously are 



10 



orthogonal. This is crucial if we want to use these modes as vectors on which we 
can project the state of the linear system PU)) and therefore obtain a simplified 
description. 

To define this scalar product consider eq. p2|) with two solutions 



E'J{x) + uj^Ej (x) + auj]Qjb{x -x^)^^ (33) 

E'l{x) + ujlEi{x) + aujfQ,d{x - Xa) = 0. (34) 

To show orthogonality the equation for £", is multiplied by Ej and vice-versa. 
Subtracting the resulting equations one obtains: 

E,E'; - EjEl' + (^2 _ u;f)E,Ej + a6{x - Xa) {E.to^Q^ - E^ufQ,) = (35) 

After integration the resulting equation on the domain [0; I] the first two terms 
drop out because of the boundary conditions. Substitute Qi, Qj in the last term 
using (PU)) leads to the the following 



(-|--n 



/ EiEjdx + ujlQiQ-j 
Jo 



0. (36) 



This shows that for toi =/= luj the term in the brackets should be zero. The scalar 
product is then defined as 

(V,; V,) = / E,E,dx + ojIQ^Q, . (37) 

Jo 

The relation ([37)) establishes a strictly positive linear form, so it is a scalar 

product. 

When the film has a finite thickness h the Dirac distribution needs to be 

replaced by the characteristic function 

g{x) — 1 for |a; — Xal < h/2, elsewhere. 

The scalar product becomes 

(V,;; V,) = / [E,E, + ulQ,Q,]g{x)dx. (38) 

Eq. p6|) shows the orthogonality of the eigenvectors Vi for the scalar product 
defined by eg. ([57)1 . Now it is possible to choose Ai such that the vectors are 
normalized 

(V,;V,) = 1. 
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Figure 3: Plot of the E component of the normahzed first, second and third 
eigenmodes. The film is in the middle of the cavity Xa = tt/2. The other 
parameters are the same as in Fig. [T] 



For this we compute (V^; V^). 

f-^ -.T ^ ,2 f' f^ ■ 2 f \ 1 sin (tOiXa) 

(V.; V.> = AU sm\u:,x)dx+ . \ '- 

\Jo sm (u;i{l-Xa)) J a 



"{uJiQ — x))dx + loq 



2-2/ N 

a sm {(jJiXa) 



~ 'V2~ 4lj, sin2Lj,(;-2;a) V2 " 2 " 4wi J {uj^ - Luf)^ 

Consequently, V^ is an orthonormal basis when 

/ X, 



. _ I Xa Sm2uJiXa 



Sin LUiXa 



Sm2uJi{l — Xa)\ LUga" Sin~ UJiXa 



Sm'^UJi{l — Xa) 



(cjg - UJ. 



2_,.,2)2 

(39) 



The normalized eigenvalues are plotted in Fig. [Altogether with the associated 
Qi in Fig. 3] for a film placed in the center of a cavity of length I = it. Notice 
the clear break in the derivative at Xa- The orthogonality of the modes (V^, Vj) 
comes from the compensation of the integral of EiEj by the product QiQj. 

These modes are specially adapted to describe the coupled system film/cavity. 
Many standard sine Fourier modes are necessary to get a good approximation of 
the first normal mode. This is seen in Fig. [5] which shows the amplitude square 
of the sine Fourier coefficients of the first normal mode. Notice the typical 1/n^ 
decay due to the fact that the second derivative of E is singular a,t x — Xa [22] ■ 

When the film is shifted to one side of the cavity, the modes become asym- 
metric as shown in Fig. [51 Again the break in the derivative is clearly apparent. 
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Figure 4: Plot of the Qi components of the normalized eigenmodes as a function 
of uj for Xa — 7r/2, 1 — tt. The eigenvalues are indicated by the vertical lines. 
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Figure 5: Plot of the squares of the Fourier sine coefficients of the E component 
of the first normal mode in log-linear scale. The function l/?i^ is plotted in 
continuous line. 
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Figure 6: Plot of the E component of the normahzed first, second and third 
eigenmodes. The film is shifted to the left of the cavity Xa = tt/A. The other 
parameters are the same as in Fig. [T] 

Here standard Fourier modes only appear for n = 4,8, ... The Qi decay very 
quickly to as shown in Fig. \7\ 

5 Cavity mode transfer using an active film 

The normal modes defined in the previous section define a basis to describe the 
state of the combined system cavity/film. We now show that it is possible by 
acting on the film to switch the cavity from one mode to another neighboring 
mode. This feature is impossible for a single linear system. It is possible here 
because of the combination of the two linear subsystems: the cavity and the 
film. 

In order to describe analytically this process, we introduce the forcing of the 
film as f{q, qt, t) and write the system as 



ett 



lit + ^oq 



~a5{x - Xa)qtt, 
ote{xa) + f{qt,t), 



(40) 
(41) 



Using the vector notation V of the previous section, this system can be 
formally written as 

{d^ + L)V = F, (42) 
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Figure 7: Plot of the Qi components of the normalized eigenmodes as a function 
of uj for Xa — 1/4^, I — TT . The eigenvalues are indicated by the vertical lines. 

where the operator L is given by (|31|1 and the forcing is 

p ^ / -aSix-Xa)f 

"V / 

For this linear system, it is natural to expand the state vector V in terms of 
the (normalized) normal modes 

■ e ■ 



(43) 



V 



E"»wv^ 



where the normal modes V^ verify the relation LV^ = ^f V^ . Plugging 
the equation (I42p and projecting over each normal mode we get 

a^ + a^Luf ^<ViF >, i = l,2,... 

where the scalar product 

<V,F>= / E,{-afSix~Xa))dx + ulQJ^iufQJ. 
Jo 

We now assume that / consists in a damping and forcing term 

/(t) = -7(i)(7, + /(i). 



(44) 



into 



(45) 



(46) 



(47) 
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Recalling the linear combination (|^H) qt = X]n cinQn we get the final expression 
of the scalar product (|1S)) 

< V,F >= -7(i)t^- a ^ AnQn + cofQJit). (48) 

To ilUustrate this, consider just two modes in the expansion (|44|) . The system 
describing the evolution of the mode amplitudes is then 

ai+aicj2 = --f{t)ujlQi{aiQi+a2Q2) + ^lQiI{t), (49) 

aa + a2Ujl ^ -7it)ujiQ2iaiQi + a2Q2) + ojlQ2l{t) (50) 

Notice that only using the normal modes (|19l20p and the scalar product (|37p 
does one obtain a consistent modal description of the system. Using for example 
the standard Sturm-Liouville modes associated to a linear impurity placed at 
X = Xa leads to an inconsistency. This important fact is shown in appendix A. 
Also remark that for I ~ t:, a — 1/2, the normal modes include the even 
(standard) sine Fourier modes. These however are decoupled from all the other 
modes because for them Q^ = so the right hand side of the amplitude equation 
(l49l) is zero. 



6 Numerical simulations 

To test these ideas, we have undertaken numerical simulations of the equations 
pn)) using the method of lines, where the spatial operator is integrated over 
reference intervals (finite volume method). The time evolution is then done 
using an ordinary differential equation solver. The algorithm is described in 
appendix B. 

6.1 Linear regime 

We introduce a characteristic forcing time function 



tanh ( I — tanh 



2 L \ wt / \ wt 

and assume that the damping and forcing are 



(51) 



l{t) - 5(i)7o, I{t) = g{t)Psm{ujt + cf>), (52) 
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Figure 8: Time evolution of the absolute difference |ai(i) — a^(t)|, i = 1,2,3 
(from left to right) between the solution of the amplitude equations (UHl) and 
the coefficients obtained by projecting the solution of the full partial differential 
equation (|10|1 . The system has been forced at frequency uj ~ 1.4 « 0^2 during 
the time 50 < t < 100. 



where ^q,(3,uj and (j) ^-re parameters. We consider the case of a centered film 
Xa — l/'2- = tt/2. For all the runs presented, we chose 70 = 0.1, /3 = 0.2,0 = 
and a time interval of forcing [ti,t2] = [50, 100] with wt ^ I- Unless otherwise 
specified, the initial state of the system is the first normal mode. 

As a first step, we validate our mode amplitude differential equations (|49p 
by comparing their solution with the mode amplitudes obtained by projecting 
the solution of the partial differential equation (fTO|) onto the normal modes V, 
given by (|19l20p . The integrals are calculated using the trapeze method using 
800 mesh points. Fig. [5]shows the absolute error in log scale as a function of time 
for the first 3 modes (except even Fourier modes). The difference is consistent 
with the error made in the trapeze integration method 0(h^). Notice that for 
1 = 2 the difference increases during the forcing. This is due to the appearance 
of new modes as shown below. The agreement is excellent and the error of about 
10~^ is essentially the error in the trapeze method h^ = (tt/SOO)^ ~ 10^^. If 
only the first three modes are used in the amplitude equations, the error is still 
very small. In many other cases, we compared the solution of the full problem 
pU]) with the one given by the amplitude equations ([i^]) and always found errors 
of about 10^^. This shows that these simple amplitude equations are a precise 
way to describe the complex system film/cavity. 

After this validation, we examine the role of the forcing frequency and show 
that we can transfer energy from one cavity mode to another by acting on the 
film via the forcing and damping (|52p. The value of the forcing frequency is 
essential as shown in the following pictures. First we chose a; — 1.4 « (^2 so 
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Figure 9: Plot of ai{t) for i = 1 (green online) and i — 2 (continuous line) after 
having forced the system at frequency lu — lA ^ lu2. The system is started in 
the mode 1 only, ai — I. 

that we are forcing the system to resonate in the second normal mode. The 
energy transfer is then efficient and after a short time of forcing we find the 
system in the normal mode 2 with very little left of the normal mode 1. This 
is shown in Fig. [9] where we plot the initial mode in dashed line and the newly 
generated mode in continuous line. This notation will be used throughout this 
section. We now change the forcing frequency to a; = 1 and retain otherwise 
the same protocol. In this case we do not have a resonance of the system and it 
responds by generating components on the neighboring normal modes. When 
forcing the system on the sine Fourier mode u^ = 2 for which Q3 = no energy 
is fed into this mode as expected from the amplitude equations (|49l) . The system 
responds on the 1st, 2nd and 4th normal modes. We have forced the system 
at frequency uj — 5 and obtained conversion from mode 1 to mode 6. This is 
shown in Fig. [TOl If we choose oj closer to oue the transfer is even better so that 
the amplitude of the mode 6 is larger. Note that it is also possible to obtain 
down conversion of modes. For example starting with mode 6 and forcing at a 
frequency uj — 0.66 « oji we obtain the mode 1 and a little of mode 2. This is 
shown in Fig. [TT] The value of Qi determines the efficiency of conversion to or 
from the mode i. For example with a cavity of length I ~ n and a film placed at 
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Figure 10: Plot ol ai{t) for i = 6 (continuous line) and i — I (dashed line) after 
having forced the system at frequency lu ^ loq — 5.06.... The system is started 
in the mode 1 only ai — 1. 
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Figure 11: Plot of ai{t) for i = 6 (dashed line) and i — 1,2 (continuous line) 
after having forced the system at frequency lo ^ lui — 0.66. The mode 3 is also 
present as shown by the blue line. The system is started in the mode 6 only 
ae = 1. 
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a = 1/2, wc find Qi = 1.8 and Qa = —0.15. We can then state that in general, 
conversion from modes close to Uq to normal modes far from ujq is more efficient 
than the converse. This is because in the amplitude equations (H^ the damping 
of Ui is proportional to Q^ while the amplification term is proportional to Qi. 
So a mode close to resonance with a large Qi is damped more than another 
normal mode with a smaller Qj. This is shown in Fig. 1101 and [TT] 

To summarize we have shown that this linear system can convert energy 
from one normal mode to another. This was thought impossible for a linear 
system because of the orthogonality of the normal modes. Here because we act 
on the sub-system we are able to do this transfer. Another important result is 
the excellent agreement between the solution of the amplitude equations and 
the solutions of the initial problem. This simple method could then be used in 
practice to solve the singular partial differential equation (fTU)) . 

6.2 Nonlinear regime 

Another way to convert energy from one mode to another is through nonlin- 
earity. A well known example is the famous study of Fermi, Pasta and Ulam 
(see for example the entry in [27j ) showing energy recurrence between Fourier 
modes in a chain of anharmonic oscillators. Now we consider the film to follow 
a law with a cubic nonlincarity and take out the driving and damping terms. 

The film equation now incorporates a cubic nonlincarity so that the compos- 
ite cavity/film is described by the system ([5]) . The cubic term can be treated as 
in the previous section and incorporated into the F term of equation (|42p . The 
scalar product is 



V,F) = LufQ.q^ = ufQ, > ^ a„Q„ . (53) 



Then the amplitude equations are 

ai+aiujl = LolQiiJ^anQnf, (54) 

n 
d2+a20jl = 0jlQ2C^anQnf, (55) 



If there is in addition forcing and damping on the system, one needs to add to 
the right hand side of (|54p the terms on the right hand side of 



20 




100 



Figure 12: Plot ofai{t) for i = 1 (continuous line), i = 2 (long dash), i = 3 (short 
dash) for a system with cubic nonlinearity. The results from the amplitude 
equations ((54l) with 8 modes are also plotted and they superpose exactly. The 
initial velocity is 0.5. The system is started in the mode 1 only. 



As the amplitude of the film polarization q increases, one expects that higher 
and lower frequency normal modes will be excited. This coupling to the other 
modes is clear from the right hand side of the amplitude equations (|M|) . Of 
course one should not increase too much the amplitude of the forcing because 
then the wave length of the cavity modes would reduce and become comparable 
with the film thickness. Then approximating the film by a Dirac distribution 
would not make sense. 

We start the system in the first normal mode with zero amplitude and a 
positive velocity and let it evolve from 0. This procedure is chosen so as not to 
create a shock in the system with the nonlinearity. For small velocities, there 
is little transfer from the first to the second and third modes. Again the sine 
Fourier modes do not play any role because for them Qi = 0. Increasing the 
initial velocity increases the transfer of energy to the higher modes. For an 
velocity of 0.5, Fig. [12] shows the evolution of a^, ?' = 1 — 4 as a function 
of time for both the full system (??) and the amplitude equations ^^. As 
expected the nonlinearity generates higher frequencies. Notice the excellent 
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Figure 13: Short time evolution of ai{t) for i = 1 — 6 from top left to bottom 
right. The same time evolutions for t — fctj-cc where tj-oc ~ 290 and k — I ~ A 
are plotted on the same panels showing recurrence. The initial velocity is 0.5 
and the system is started in the mode 1. 

agreement between the solution of the partial differential equation system and 
the amplitude equations. This holds even for such large velocities as 2. 

We conclude this section with an observation of recurrence similar to what 
happens for the Fermi-Pasta-Ulam system. Fig. [13] presents a short time evo- 
lution of the first 6 modes for an initial velocity of 0.5, starting from the first 
mode. In the plots of Fig. [13] are superposed 4 other time evolutions taken at 
times t = fcircc where irec — 290 and k is an integer. This recurrence could 
indicate that our system is close to being integrable. 

7 Conclusion 

We considered the interaction between an electromagnetic field in a cavity and 
a thin polarized dielectric film. The model is the Maxwell-Lorenz system where 
the medium is described by an oscillator and the coupling to the wave equation 
is through a Dirac delta function. We introduced normal modes which are 
adapted to the system film/cavity. These are well adapted to describe the 
time evolution of the system, unlike the standard sine Fourier modes or other 
Sturm Liouville modes. The normal modes composed of the field E and the 
displacement Q are orthogonal with respect to a special scalar product which 
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we introduce. The amplitude equations derived from the normal modes give an 
excellent description of the dynamics and could even be used as a numerical tool 
instead of solving the full partial differential equation system using the fairly 
involved finite volume method. 

Assuming a linear oscillator for the film, we show conversion from one normal 
mode to another by forcing the film at specific frequencies. This is new for linear 
systems and could be used for many applications in optics or microwaves. 

If the film is described by an anharmonic oscillator, the evolution generates 
other modes. Again the amplitude equations provide excellent agreement with 
the solution of the full problem. Finally we observed recurrence for certain 
initial velocities of the film. This phenomenon is known to exist for systems 
close to integrability. The fact that we observe it here may indicate that our 
system is in some ways close to integrability. 
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7.1 Appendix A: Inconsistent projection using standard 
eigenmodes 

Using standard Sturm-Liouville eigenmodes and the usual scalar product leads 
to inconsistent results. To show this let us assume no forcing for simplicity. We 
consider the usual eigenmodes associated with the Sturm Liouville problem 

E^^-a'^E6{x-Xa)^-uj^E. (56) 

Call these modes En associated to the eigenfrequency cun- 
One then expands the field as 

oo 

e(x,t) = y^^an{t)En(x). 

Plugging this expansion into the system of equations (jlOp and projecting onto 
the En using the standard scalar product, one gets the evolution of the a's 

di + Lo\a\ — aujlqEi{a), (57) 

d2 + a;|a2 = aujlqE2{a), (58) 

q + ujf^q = a^aiEi{a). (59) 
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Let us examine the jump condition on E^ from the partial differential equation 
CT- We have 

- [E^TJ- = -ai^uj^q + aE{a)]. (60) 

From the expansion (17. ip . we obtain 

-[-B^r]""^ = -S^ aniEnJ'^l. 
n 

We get the jumps of the Enj. from the eigenvalue relation 

Enxx ~ a^En5{x - Xa) = -Uj'^En, 

and obtain that 

-'^ciniEnxTj'- = -a'^'^a^Enia) = -a^E{a), 

n n 

which is clearly inconsistent with the result obtained from the original system 
(pn|) . Therefore one needs to use the normal modes associated with the full 
system and the special scalar product (p7|) to get a consistent reduced description 
of the system. 

7.2 Appendix B: Numerical method for solving ( 1101) 



The evolution equation pop involves a partial differential equation with a Dirac 
distribution coupled to an ordinary differential equation. We solve this coupled 
system using the method of lines where the time operator is kept as such and the 
space operator is discretized, naturally leading to a system of coupled ordinary 
differential equations. The spatial operator is a distribution so the natural 
way to give it meaning is to integrate it over small reference intervals (finite 
volume approximation). The value of the function is assumed to be constant 
in each volume. This method of lines allows to increase the precision of the 
approximation in time and space independantly 

We tranform pI7|) into a system of first order partial differential equations 

et = /, (61) 

ft=exx-aS{x-Xa)rt, (62) 

qt = r, (63) 

rt = -Lolq + ae(xa). (64) 
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We then define our volume elements making sure that Xa is the center of one 
element. We then integrate the operator on each volume [a;„ — ft,/2,x„ + h/2\ 
where h is the space step. For Xn ^ Xa, we recover the standard finite difference 
expression for the Laplacian 

Cri+i + e„+i — 2e„ 



^2 

and get the discrete wave equation 



0{h'), 



e„ = /„ (65) 

Cn+i + e„+i — 2e„ 



/l2 



(66) 



For Xn = Xa: we get 



e„ = /„ (67) 



/l2 



ar (68) 

= r, (69) 



f ^ -uj^q + aen- (70) 

The coupled system (|65l67p is then integrated using an ordinary differential 
equation solver. In practice we use the variable step Runge-Kutta 4-5 Dopri5 
developped by Hairer and Norsett at the University of Geneva [^SJ . 
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